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Abstract. With the understanding that the enigmatic Gamma-Ray Burts (GRBs) are beamed explosions, and with the recently 
discovered "Ghiiianda-relation", the dream of using GRBs as cosmological yardsticks may have come a few steps closer to 
reality. Assuming the Ghiiianda-relation is real, we have investigated possible constraints on cosmological parameters using a 
simulated future sample of a large number of GRBs inspired by the ongoing SWIFT mission. Comparing with constraints from 
a future sample of Type la supernovae, we find that GRBs are not efficient in constraining the amount of dark energy or its 
equation of state. The main reason for this is that very few bursts are available at low redshifts. 
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1. Introduction 

The usage of thermonuclear Type la supernovae (SNe la) 
has revolutionized cosmology. These intrinsically bright ex- 
plosions are almost standard candles in optical light. With 
a simple light-curve correction they can be standardized to 
high enough precision to probe in detail the energy content 
of the universe. This led to the discovery that dark en ergy 
dominates the presently acce lerating universe (e.g. jRiess et alJ 
1998; Perlmutter et al. 1999). Today, SN la data is approaching 
the quantity and quality where it is possible to constrain not 
only the amount of dark energy but also dark energy properties 
(e.g. JHannestad & Mortseli2002lEo04l) . 

Gamma-Ray Bursts are even more powerful explosions. 
They have recently been firmly linked to energetic core - 
collapse supernovae ( H iorth et al. 2003: Mathes on et al.l2 003). 
Their isotropic energy appears to outpower thermonuclear su- 
pernovae, which allows studies at even higher redshifts. Also, 
in contra st to the case of SN e la where the rate is unknown at 
z > 1.5 (Dah len et all2004. we know that GRBs exist at high 
redshifts ("e.g.. lAndersenet aljEo OO). Moreover, the released 
burst of gamma-rays can penetrate the dust that obscures our 
view of the distant universe in optical lig ht. For current con- 
straint s on dust attenuation of SNe la, see Ost man & Mortseill 
J200.4 . The potential of GRBs as probes for cosmological in- 
vestigations thus appears to be very good. 

The realization that the total gamma-ray energy of a GRB, 
when corrected for the e ffects of beaming , spans a reasonably 
narrow range of energies ( Frail et al. 2001 ) arose hope for GRB 
cosmology (e.g.. lSchaefeJ2003l) . More recently, the discovered 
tight relation between the rest-frame peak energy E pelL k of the 
GRB and the rest-frame, beaming-corrected ga mma-ray energy 
release E y , the so-called "Ghirlanda-relation" jGhirlanda et alJ 



2004a), has renewed the hope for GRB cosmology. This rela- 
tion allows an empirical correction to the determined luminos- 
ity distances for each GRB, in much the same way as light- 
curve shape corrections are applied to SNe l a, with a scatter 
along t he GRB Hubble diagram of ~ 0.5 mag ( Ghirl anda et alJ 
l2004bl) . This has caused a f ierce activity of research in the 
area of GRB cosmology (e.g..lGhirlanda et alJ2004bllDai et all 
2004; iFirmani et all2005t IXu et all2005HXul2005h . 

There is, however, much to be done before GRB cosmology 
can be established. The reali ty of the Ghi r landa relation is still 
under discussion dBand & Preecd l2005t IFriedman & BloomI 
2005t iGhirlanda et al.l l2005). Moreover, very different results 
on GRB cosmology are obtained by the different authors, based 
on the current sample of GRBs. Friedman & Bloom] d2005l) 
thoroughly examined pro and cons of GRB cosmology with 
the presently limited sample of well studied GRBs, and show 
that the obtained results are crucially dependent on the choice 
of (poorly known) input parameters. The way to select w hich 
GRBs to include also influence the results. IFriedman & BloomI 
d2005l) are therefore rightfully cautious concerning claims of 
the utility of GRBs for cosmology. A larger dataset is clearly 
needed to establish this issue. 
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The recently launched SWIFT satellite JGehrels et alJ2 004) 
will find hundreds of GRBs. In this paper we simulate the 
potential effect of such a large number of GRBs for the use 
of determining cosmological parameters, assuming that the 
Ghirlanda relation holds. 

In Sect. |2 we present the method used for deriving cos- 
mological parameters from GRBs. We also discuss various as- 
sumptions for the constructed input samples that we use for the 
simulations. The results are discussed in Sect.|3] 



2 



E. Mortsell and J. Sollerman: GRB cosmology 



2. Simulations 

2.1. The method 

In the follo wing, we have basically fol lowed the formalism laid 
out in, e.g.. lFriedman & Blooml J2005). GRBs are used in much 
the same way as SN la standard candles in constraining cos- 
mological parameters, i.e, by comparing observed luminosity 
distances to theoretical predictions. However, since the process 
of standardizing the GRB candles is cosmology dependent, the 
Ghirlanda relation needs to be recalibrated for each cosmol- 
ogy. Effectively, this amounts to refitting for each cosmology 
the £peak - E y power-law 

IE V 

£peak=W-^l ■ (1) 

Note that Eq is an arbitrary constant and that by putting 
Eq oc /i~ 3 / 2 , the best-fit k and rj will be independent of the 
value of the Hubble parameter. Thus, marginalizing over k 
and T] is similar to margin alizing over M for SNe la (e.g., 
IHannestad & Mortsell2002l) . 

There have been several suggestions on how to include 
the informati on fro m the £ P eak - E Y fit in the cosmology 
fit. iDai et alJ {2004) basically ignore d this complicat i on. Thi s 
was quickly noted and remedied bv iGhirlanda et al ] d2004bl), 
who ref itte d the relation for each cosmology. iGhMandaetall 
( 2004bjl and lXu et alJ (12005 ) have also used a simple, but rather 
unrealistic, approach in assuming that k and 77 can been exactly 
determined using a sample of low redshift GRBs, or by the- 
oretical co nsideratio ns, an d does not need to be recalibrated. 
iFriedman & Blooml J2005b also refits the £p ea k - E y power-law 
and obtains a new set of k and 77 with corresponding errors for 
each cosmology. These errors are then propagated to the er- 
ror in the derived luminosity distance. This approach has the 
drawback of giving smaller ^2-values for cosmologies where 
the £peak - E y re lation is ba dly fit since the luminosity distance 
errors are larger. Xu (2005) discuss a method where cosmolo- 
gies with a good fit of the Epeak - E y relation are favored (their 
Method III). This method is similar to adding the xl- values 
from the cosmology fit and the power-law fit for each cosmol- 
ogy- 

In this paper, we treat k and 77 as unknown parameters that 
should be marginalized over. This means that we do not use 
the k and 77 that gives the best fit to the £p ea k - E y relation for 
each cosmology, but instead use the k and 77 that minimizes the 
cosmology ^2-value, in analogy with marginalizing over Al for 
SNe la. We have noted that our constraints on the cosmological 
parameters does not depend sensitively on the specific method 
used, as long as k and 77 are not assumed to be fixed. 

2.2. The GRB sample 

The aim of this paper is to investigate the future potential of 
GRB cosmology by simulating a larger sample of GRBs. For 
supernova cosmology, many such investigations have been per- 
formed, i n particular in conjunction with the plan ned SNAP mis- 
sion (e.g- lGoiiath et al.l200lllWang et al.l2004 . 



To simulate the effect of a larger future sample we have 
chosen 200 GRBs with essentially similar properties as the cur- 
rent sample. This number of GRBs is in line with the expecta- 
tions on the SWIFT satellite, which is predicted to find about 
100 b ursts every year for a life time of 2-8 years JGehrels et alJ 
2004). 

2.2.1 . The redshift distribution 

An important property of the future sample is the redshift dis- 
tribution of the GRBs. It has been suggested that GRBs and 
their afterglows can be detected up to very high redshifts ( Lamb 
I2002I) . To investigate the dependence of the cosmological pre- 
dictions on the distribution of the objects we have used two 
published predictions for the redshift distribu tion. 

First, the calculations performed by Bromm & Loeb (2002) 
suggest that the SWIFT satellite will be able to detect GRBs up 
to redshifts of z > 20. They find that 25% of the detected bursts 
will have z > 5. We have adopted their calculated redshift dis- 
tribution for SWIFT as input for our simulations. 

Secondly, a more conservative estimate is performed by 
iGorosabel et alJ yP04). They argue that relatively few high- 
redshift GRBs will be found by SWIFT, and calculating the ex- 
pected redshift distribution from their Fig. 3 we distribute the 
simulated GRBs up to z ~ 6. These two distributions, as shown 
in Fig. n can thus be taken to represent two extremes for a 
future GRB sample. 

We note that systematic effects such as gravitat ional lens- 
ing a r e potentially greater at high redshifts feergstrom et alJ 
2000; lAmanullah et aDl2003h . In this study, we assume these 
errors to be negligible. Note also that the probability for mul- 
tiple lensing increases with redshift and that multiple im- 
aged GRBs therefore are potentially useful for consteaining 
the Hubble parameter and galaxy halo properties jGoobar et alJ 
200i lMortsell et all2005l) . 

2.2.2. The properties of the sample 

We have taken the current sample of 19 GRBs from the list pro- 
vided by www.comsicbooms.net, and distributed the observ- 
ables and their associated errors in the simulated sample in the 
same way as for the observed sample. 

The fractional errors are, in order of decreasing importance 
for the total magnitude error, the peak energy 0" £ /£p e ak = 
0.176, the jet break time cr,/t = 0.18, the cirumburst medium 
density cr„/n = 0.5, the fluence os/S = 0.1 and the k- 
correction cr^/k = 0.056. The contributions to the magnitude 
error are 0%°* ~ 0.4, cr' m ~ crl ~ 0.2, cr s m ~ 0.1 and 
oj n ~ 0.06. The total magnitude error is cr m ~ 0.5. 

We have chosen to assume that the uncertainties in the sim- 
ulated future sample will be the same as in t he current sampl e 
of GRBs. This is similar to the approach bv lXu etafl J2005). 
We also assume that even the most distant bursts will be fol- 
lowed in enough detail not to deteriorate the sample. In the 
gamma-ray regime, it is clear that the intrinsic scatter in pa- 
rameters is already larger than expected from the redshift dis- 
tribution alone. 
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For a true SWIFT sample we should distribute the £p e ak m a 
smaller range, since the BAT instrument is sensitive only in the 
relatively narrow 15-150 keV range. We have, for the sake of 
simplicity, ignored this matter, but note that it has been used as 
an argument to continue the efforts for HETE II a nd Integral 
jFriedman & BloorrJ2005llGorosabel et all2004l) . 

iLambl J2002I) argued that even for the most distant GRBs, 
follow-up observations will be possible. Redshift determina- 
tions will indeed be feasible up to z ~ 10 with instruments such 
as X-shooter (D'Odorico et al. |2004 . While optical follow-up 
of very distant bursts may not be good enough to probe the jet 
break in detai l, also X-ray and near-IR facilities are available 
for this. iGorosabel et alJ J200 4) estimate that a near-IR after- 
glow can be followed up to z = 9 even with modest exposure 
times on a 10-m class telescope. The very late occurrences of 
jet-breaks will of course be difficult to detect for the dimmest 
targets, which may bias the sample. 

The narrow energy range for SWIFT means that 200 GRBs 
with measurements of all relevant observables is probably a too 
optimsitic assumption for this mission. This is also supported 
by the fact that the SWIFT bursts detected so far appears to 
be faint and has the refore not been successfully characterised 
tBerger et alJl2Q05l) . The sample we use for our simulations 
should therefore be regarded as a rather optimistic guess that 
has been influenced by the SWIFT mission, but may have to 
await future missions to be accomplished. This optimsitic as- 
sumption will only strengthen our conclusions given below. 

2.3. The SN sample 

To compare the constraints on the cosmological parameters ob- 
tained from our GRB simulation we have also simulated a fu- 
ture set of SN la data usi ng the publicly available SNOC pack- 
age (Goo bar et alJ l2002). Just as the constructed GRB sample 
was motivated by the SWIFT mission, which is already ongoing, 
we base the constru cted SN la sample on the available Gold 
sample jRiess et al.ll2004l) as well as two ongoing supernova 
surveys, the ESSENCE project and the Supernova Factory. 

The ESSENCE project (e.g.. iMatheson et allEooll) is an 
ongoing survey aimed to measure 200 SNe la in the redshift 
domain z = [0.2 - 0.8]. The goal is to derive tight constraints 
on the equation of state, w, of the dark energy. In construct- 
ing the sample for our simulations we have used the redshift 
distribution of the hitherto discovered 109 SNe la, but in- 
creased the number of SN la to 200, as planned for ESSENCE 
(Mik naitis et al J 120041) . We have adopted an intrinsic photo- 
metric error for an individual supernova of 0.14 mag. Binning 
200 SNe into Am = 0.1 redshift bins gives a statis tical un- 
certainty JGarnavich et al.l2002HMiknaitis et al.l2005l) close to 
the systematic floor expected from e.g., uncertainties in the k- 
corrections. 

The Supernova factory is an ongoin g effort to measure 30 
nearby, z = [0.03 - 0.08], SNe la (e.g.. lAldering et alJl2002t) . 
We have assumed an uniform redshift distribution in this range, 
and again an individual error of 0. 14 mag per supernova. Note 
that our approach only includes statistical errors. 



3. Results and discussion 

In Fig. 13 we provide the results of our simulations. All simu- 
lations have been made assuming a flat universe (Qm = 0.3, 
£2a = 0.7) dominated by a cosmological constant (w = -1). 

In Fig. |2 we show the probability contours for the 
200 GRBs (dashed lines) and the combined (Gold sam- 
ple+ESSENCE+SN factory) SN la sample (solid lines). 
Combined results are showed using yellow (filled contours). 
The four contours correspond to 68.3 %, 90 %, 95 % and 99 % 
confidence, respectively. 

From the top row in Fig. |2 it can be seen that the GRB 
sample (dashed lines) in itself is not very efficient in constrain- 
ing Qa. This is true for both adopted redshift distributions. The 
interesting aspect is that the constraints are rather complemen- 
tary to the SN la constraints, which makes the combination of 
these two datasets relatively fruitful. That the GRB cosmology 
is predominantly sensitive to Am, and rather insensitive to the 
value of the cosmological constant, reflects the higher redshift 
distribution of this sample. At these epochs, the dark energy 
contribution to the energy density was still relatively unimpor- 
tant and we lack a set of low-z GRBs that would provide the 
necessary leverage in the Hubble diagram to constrain the dark 
energy. 

Clearly, a good prior on Q M from large scal e structure 
LSS) surveys such as 2dF JPercival et alJl200lh or SDSS 
Tegmark et al. 2004) are more helpful in constraining the SN 
la contours, but we find it noteworthy that the same thing can 
be achieved using only standard candle techniques. 

The lower row in Fig. |2] shows constraints derived for a 
constant equation of state parameter, w = wo, for the dark en- 
ergy. When fitting wq, we have assumed a flat universe, i.e., 
Hm + Qx = 1 ■ The solid contours correspond to the SN con- 
straints and the dashed contours to GRB constraints. From 
this exercise it can be seen that GRB cosmology can not be 
expected to give very useful constraints on the equation of 
state p arameter. It has been suggested (e.g., Frie dman & Blooml 
2005) that the GRBs could be useful to constrain any potential 
evolution of w. We believe that, given the large uncertainties 
in even constraining a constant value for the equation of state 
parameter, such efforts will be in vain. It is interesting to note 
that combined with a good prior on Qm from, e.g., LSS sur- 
veys, SN la data alone will provide very powerful constraints 
on a constant equation of state parameter within a few years. 

Our main conclusion is thus that GRBs are not likely to 
contribute significantly to constraints on the cosmological 
paramete rs in the near f uture. This is similar to the wordings 
in Friedman & Bloom ( 2005) and the also to the results 
of IXu et alJ IJ2005) when not assuming an artificially fixed 
£p ea k — E v relation but at odds with other investigations 
(e.g., Ghis ellini etal] 120051: lLazzati et all l2005h . The main 
reasons for the limited use of GRBs, apart from the huge 
inher ent uncertainties in trying to standardize these candles 
(e.g., iFriedman & Bloomll2005t) . is their redshift distribution. 
Probing Ox and its equation of state benefits from a wide 
redshift distribution, with good coverage at the epochs where 
£2x dominates the energy density and a set of low-z events 
to constrain the normalization of the GRB luminosity. By 
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artificially including a large number of low-z bursts into 
the sample, we were indeed able to considerably shrink the 
confidence regions for the cosmological parameters. The 
point is that the calculations used for the future GRB redshift 
distributions used in this work predict that very few of these 
rare bursts will be available within the limited volume spanned 
by these low redshifts. 
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Fig. 1. The redshift distributions used for our simulated GRB 
sample. The dashed line is the predicted redshift distribution 
available to SWIFT according to lGorosabel et alJ d2004l) . while 
the solid line s hows the distribution of GRBs according to the 
calculations by Bro mm & LoebH2002l) . 
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Fig. 2. Upper left: Q . \ vers us Qm for the GRB sample based on 
the iGorosabel et all 1 20041) redshift distribution (dashed lines) 
and the SN sample (solid lines). Combined results are showed 
in yellow (filled contours). Upper right: Qa versus Qm for 
the GRB sample based on the Brom m & Loebl J2002I) redshift 
distribution. Lower left: wq versus Q M assuming a flat uni- 
verse for the GRB sample based on the Gorosa bel et alJ (|2004) 
redshift distribution (dashed lines) and the SN sample (solid 
lines). Combined results are showed in yellow (filled contours). 
Lower right: wn versus Q M for the GRB sample based on the 
Bromm & Loe hi J2002h redshift distribution. 



